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Variability in the Lambda Orionis cluster substellar domain* 
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Abstract. We present the first results on variability of very low mass stars and brown dwarfs belonging to the ^5 Myr 
Lambda Orionis cluster (Collinder 69). We have monitored almost continuously in the J filter a small area of the cluster 
which includes 12 possible members of the cluster during one night. Some members have turned to be short-term variable. 
One of them, LOril67, has a mass close to the planetary mass limit and its variability might be due to instabilities produced 
by the deuterium burning, although other mechanism cannot be ruled out. 
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1. Introduction 

We are embarked in a very ambitious project devoted to the 
Lambda Orionis Star Forming Region (SFR), the Head of the 
Orion giant, the Hellenic mythological hero. This area con- 
tains several distinct structures, including several dark clouds 
with active star formation and a very young cluster (Collinder 
69 or Lambda Orionis cluster) associated to the multiple star 
A Ori, located at the center of the SFR. Additional details can 
be found in Dolan & Mathieu (1999, 2001) and Barrado y 
Navascues et al. (2004, 2005). This paper deals with a spe- 
cific topic: variability at the bottom of the stellar sequence in 
the Lambda Orionis cluster. Specifically, we are interested on 
the faintest object in our sample, trying to establish whether 
the variability arises from pulsation due to deuterium burn- 
ing, a mechanism initially proposed Toma (1972) -based on 
previous work by by Gabriel (1964)- for pre-main sequence 
stars (2.0-0.2 M0) and recently suggested by Palla & Baraffe 
(2005) for young brown dwarfs. 

So far, few young brown dwarfs have been photomet- 
rically monitored, and their periods derived (see Scholz et 
al. 2005, this volume). The pioneering work dealt with a 
handful of very low mass members of the Pleiades cluster 
(Terndrup et al. 1999; see also Scholz & Eisloffel 2004a). 
The youngest objects belong to Chamaeleon I (l-3Myr), 
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and they have the longest periods, in the range of 53-114 
h and with amplitudes in the J band of 0.13-0.05 mag (Jo- 
ergens et al. 2003). Several works have focused on brown 
dwars of the somewhat older (3-5 Myr) a Orionis cluster 
(Bailer- Jones & Mundt 2001; Caballero et al. 2004; Scholz 
& Eisloffel 2004b). They have unveiled periods in the range 
3-240 h, with amplitudes of 0.02-0.4 mag. However, the low- 
est mass object whose variability has been detected so far, 
S Ori 45, has a mass about 0.02 Mq. Its period has a been es- 
timated as 0.50±0.12 h (Bailer- Jones & Mundt 2001). Later 
on, Zapatero-Osorio et al. (2003) derived two possible peri- 
ods, either ^Q.ll h or a value in the range 2.5-3.6h, with 
amplitudes in the range 0.01-0.13 mag. Because of escape 
velocity arguments, the shortest period cannot be the rota- 
tional period. Palla & Baraffe (2005) suggested the possibil- 
ity of pulsations due to deuterium burning. However, except 
for S Ori 45, and taken into account the age, the other objects 
seem to be too massive (larger than 0.05 M©) to be undergo- 
ing deuterium burning (ie., they should lack pulsations driven 
by its combustion). We will discuss the observations and data 
analysis in the next sections. 

2. Analysis 

2.1. The observations 

The data presented in this paper were collected during two 
campaigns in November 19-20 2002 and February 12-14 
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2.2 Analysis of the variability 




Fig, 1. Deep INGRID image in the J band of the area south- 
east of the star A Orionis. We have marked the cluster candi- 
dates from Dolan & Mathieu (1999) and Barrado y Navascues 
et al. (2004). 



2003 at the 4.2m WHT, La Palma (Spain), and the near in- 
frared camera INGRID, which has a 4.2x4.2 arcmin^ field 
of view. During the first run, we divided the area around the 
Lambda Ori cluster in a 42x28 arcmin^ grid, composed of 
squares of 4 x4 arcmin each, trying to match our CFHT opti- 
cal survey (Barrado y Navascues et al. 2004). We used the J 
filter, but in few cases we took additional data with the filters 
H and Ks. The exposure times were five individual expo- 
sures (with a small offset) for each pointing of 60 seconds 
each. Therefore, the total exposure time was five minutes per 
pointing. The same procedure was repeated during the 2003 
run. These observations will be reported elsewhere. The data 
were processed using IRAF' . 

Here, we analyze the data taken just southeast of the A 
Orionis star. In 2002 we collected JHKs imaging of this 
pointing, whereas in 2003 we devoted a whole night, with 
a five-hour time-span, in order to observe this small area in 
the J band, looking for short term variability. The individ- 
ual 60 seconds exposures were processed, aligned and com- 
bined into one single 5 minutes frame. In total, we collected 
51 of such five-minute images. In addition, we combined all 
the frames into a single image in order to derive deep, ac- 
curate photometry. The final deep image can be seen in Fig- 
ure 1, where we have marked the location of our targets. The 
star BD +09 879C, the third component of the A Ori system 
(30 arcsec south of the primary), can be partially seen at the 
upper-right corner of the finding chart. Optical and near in- 
frared photometry, as well as a mass estimate based on a 5 
Myr isochrones by Baraffe et al. (1998) are listed in Table 1. 

' IRAF is distributed by National Optical Astronomy Observato- 
ries, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under contract to the National Science 
Foundation, USA 
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Fig. 2. Light-curves for two low mass stars belonging to the 
Lambda Orionis cluster These objects have been classified 
as Classical TTauri stars based on their location on infrared 
color-color diagrams, using Spitzer data. The data for the 
comparison stars are shown at the bottom of each panel. 

Note that we have previously collected low resolution 
spectroscopy for two of the objects listed in Table 1, LOri 
134 and 135. The spectral types are M5 and M7, and they 
were classified as probable non-member and possible mem- 
ber, respectively, based on their spectral type and the optical 
and 2MASS photometry (less accurate than our new values). 
Based on the new photometry (including Spitzer/IRAC data, 
see Barrado y Navascues et al. 2006, in preparation), we have 
classified them as cluster members. The other objects ana- 
lyzed here, except for LOril65, seem to be members of the 
association too. 

2.2. Analysis of tlie variability 

The J-band light-curves of the bona-fide members are dis- 
played in Figures 2-A. The first diagram displays the curves 
for two low mass stars, whereas the second shows the data for 
two probable brown dwarfs. The last figure (Figure 4) corre- 
sponds to an object with a mass around the planetary mass 
limit (^14 Mjup). These light-curves have been constructed 
subtracting the magnitudes of the targets from a non-variable, 
artificial comparison star 

In order to construct the artificial comparison or refer- 
ence star we have analyzed the light-curves of field stars with 
magnitudes as similar as possible to those of the targets. From 
these objects, we have selected those that do not display sig- 
nificant variability over the five hours of observations The 
reference star is built by selecting field stars in a range of ±2 
mag around the problem star. After several iterations to re- 
move variables in the field, we retain between 5 and 7 stars 
with no apparent variability. On average, the standard devia- 
tion of their brightness is below 0.008 magnitudes. We have 
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Table 1. Photomotry of our candidate members. 



Name 


Rc error 


Ic error 


J ^ error 


error 


//■^ error 


Ks^enoi 


Mass'' (Mq) 


Sp.T. 


Var? 


[DM99] 35^ 


14.48 0.01 


13.61 0.01 


12.136 0.001 


12.26 0.01 


11.60 0.01 


11.55 0.01 


0.785 0.739 0.571 






LOriOOe 


13.55 0.01 


12.75 0.01 


11.430 0.001 


11.67 0.01 


10.90 0.01 


10.94 0.01 


1.106 1.171 0.901 






LOriOie 


14.07 0.01 


13.18 0.01 


11.902 0.001 


12.00 0.01 


11.41 0.01 


11.27 0.01 


0.970 0.875 0.696 






LOri030 


14.95 0.01 


13.74 0.01 


12.314 0.001 


12.44 0.01 


11.81 0.01 


11.64 0.01 


0.735 0.655 0.537 






LOri048 


15.78 0.01 


14.41 0.01 


12.622 0.002 


12.78 0.01 


12.17 0.01 


12.00 0.01 


0.498 0.534 0.418 






LOri063* 


16.80 0.01 


15.34 0.01 


13.610 0.002 


13.72 0.01 


13.02 0.01 


12.69 0.01 


0.287 0.269 0.261 


M4.5 


Y 


LOri085* 


17.65 0.01 


16.04 0.01 


14.058 0.003 


14.21 0.01 


13.58 0.01 


13.26 0.01 


0.184 0.194 0.170 




N 


LOril02 


18.24 0.01 


16.50 0.01 


14.430 0.004 


14.57 0.01 


14.05 0.01 


13.78 0.01 


0.137 0.147 0.113 






LOril34 


19.91 0.01 


17.90 0.01 


15.609 0.006 


15.72 0.01 


15.17 0.01 


14.82 0.01 


0.063 0.059 0.046 


M5 


Y 


LOril35 


19.91 0.01 


17.90 0.01 


15.529 0.006 


15.63 0.01 


15.14 0.01 


14.79 0.01 


0.063 0.063 0.048 


M7 


Y 


LOril65^ 


23.12 0.22 


20.73 0.02 


18.497 0.026 


18.57 0.05 


18.08 0.05 


17.83 0.09 








LOril67 


23.86 0.64 


20.90 0.02 


17.776 0.017 


17.84 0.03 


17.15 0.03 


16.62 0.03 


0.018 0.014 0.015 




Y 


^ INGRIDFeb. 12, 2003. 



2 INGRID Nov. 19, 2002. 

^ From Dolan & Mathieu (1999), the spectral type is in the range K8-M2, from the colors. 
* Class II object, based on Spitzer/IRAC data (Barrado y Navascues et al. 2006). 
^ Possible non-member. 

® As derived from /, J and K magnitudes, (m — Af)o=8.010, E{B — V)=0.12 mag and 5 Myr isochrones from the Lyon group. 



weighted them according to their observational errors and 
computed their average magnitude. We have used a Scr as 
our variability criterion. In all figures, the light-curves of the 
comparison stars together with their mean J-band magnitudes 
are shown below the light-curves of the targets. 



2,3. Variability on the context of the evolutionary status 

Figure 5 displays the HR diagram for our targets. Bolo- 
metric luminosities were derived after the / magnitudes the 
bolometric corrections by Comeron et al. (2000), a distance 
modulus of 8.010 mag and A/=0.223. Similar values can 
be derived from J and K using other bolometric correc- 
tions. Effective temperatures come from Bessell et al. (1991) 
and the (i? — /) color Note that the Tcff scale change very 
much the Premain-sequence locus. For instance, the scale by 
Leggett (1992) would move the cluster sequence close to the 
1 Myr isochrone. Moreover, error-bars for a given temper- 
ature scale can be close to 150-200 K. For the dataset rep- 
resented in the diagram, the 5 Myr Baraffe et al. (1998) fits 
the location within the errors. The figure also includes the D- 
instability strip (delimited by thick, dashed lines), which ap- 
pears due to deuterium burning during the first million years 
of pre(sub)stellar evolution. 

LOri085 and LOri063 can be classified Class II objects 
based on Spitzer/IRAC data (Barrado y Navascues et al. 
2006). For the first star, it seems that the variability, if present, 
is very small, at the same level as in the case of the com- 
parison star (in our dataset, crLOri085 =0.007 mag, whereas 
fcomp=0.004 mag). However, there is an obvious change in 
the case of LOri063. Its variability can be due either to spots 
or changes in the accretion rate. In the first cases the vari- 
ability would be essentially periodic (with a range from few 
to hundred of hours), whereas in the second case the light- 
curve would lack periodicity. Our dataset do not allow us to 
distinguish between these possibilities. 
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Fig. 3. Light-curves for two brown dwarf candidates. 



LOril34 and LOril35 seem to be bona-fide brown dwarfs 
based on photometry and their spectral type. Based on 
Spitzer/IRAC data, they can be classified as Class II and Class 
III objects, respectively, although LOril35 might have IR ex- 
cess at 3.6 and 4.5 micron. At the timescale of our photomet- 
ric monitoring, the light-curves seem to be very different. 

LOril67 is a cluster member based on its optical and 
near IR data and its spectral energy distribution (from the 
R band to 8.0 micron, Barrado y Navascues et al. 2006). It 
is located close or within the deuterium string instability in 
a HR diagram (Figure 5). The question is: Are we detect- 
ing pulsation? The light-curve suggests a period of about 
one hour. However, error-bars are large and the time cov- 
erage is not optimal to confirm this suggestive possibility. 
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Fig. 4. Light-curve for LOril67. If member, this object 
would have a mass slightly above the the deuterium burning 
limit, located at 0.013 M0. Due to the cluster age, pulsations 
might be present. 



Note that the theoretical prediction by Palla & Baraffe (2005) 
gives a period P(D-instability)^1.5 h. Moreover, for an age 
of 5 Myr, the radius of LOril67 would be 0.225-0.143 R© 
for a 0.02 Mq brown dwarf (NextGen or COND models by 
the Lyon group). The escape velocity for this object would 
be 173 or 231 Km/s for each model, after applying the for- 
mulae Voscapc=617.88x[M(M0)/R(R0)]O-5. Therefore, the 
minimum rotational period would be Prot,min=27rR/Vcscapo, 
which translates into 1.79 h or 0.75 h for for NextGen and 
COND models, respectively. Since for the mass and Toff of 
LOril67, the COND models are better suited to describe its 
properties, and due to its location in the HR diagram (very 
close or within the D-instability strip), it seems that it is a 
real possibility that we have discovered a D-burning pulsa- 
tion. In any case, as stated before, this is just a possibility, 
and other explanations, more conventional, are still possible 
(such as variability due to spots or clouds in the atmosphere 
of the object). Additional data, covering several nights and 
with better phtometric accuracy, are needed. 



3. Conclusions 

By photometrically monitoring in the J band a sample of low 
mass members of the Lambda Orionis cluster, both stars and 
substellar objects, and using differential photometry, we have 
derived the light-curves in a five-hour time-span. We have 
detected variability in some of them, including a Classical 
TTauri star, two brown dwarfs and a very faint object with 
J=17.8 mag, which corresponds to a mass close to the plane- 
tary mass domain for cluster members. This variability, with 
an amplitude of about 0.1 mag, might be due either to ro- 
tation, accretion or to instabilities -ie., pulsation- due to the 



Fig. 5. HR diagram for our cluster candidates. We have in- 
cluded isochrones (from top to bottom 1, 3, 5, 10, 20, 50, 100 
and 1000 Myr, solid lines) and evolutionary tracks (from left 
to right 0.02, 0.03, 0.04, 0.06, 0.1, 0.2, 0.4, 0.6, 0.8, 1.0 Mq, 
dotted lines) from Baraffe et al. (1998). The 5 Myr isochrone 
has been highlighted. The area delimited by the two thick 
dashed lines correspond to the deuterium instability, where 
pulsation might occur 



deuterium burning, as proposed by Palla & Baraffe (2005) for 
very low mass objects. 
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